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Abstract 

Using two different methods, we have determined the rescaling of the scalar 
condensate Z = Z v near the critical line of a 4D Ising model. Our lattice data, in 
agreement with previous numerical indications, support the behavior Z^ ~ ln(A), 
A being the ultraviolet cutoff. This result is predicted in an alternative description 
of symmetry breaking where there are no upper bounds on the Higgs boson mass 
from 'triviality'. 



1 Introduction 

There are many computational and analytical evidences pointing towards the 'triviality' 
of $ 4 theories in 3 + 1 dimensions [1,2] (see also [3] and references therein), though a 
rigorous proof is still lacking. Nevertheless these theories continue to be useful and play 
an important role for unified model of electroweak interactions. 

The conventional interpretation of these theories extends to any number N of scalar 
field components and, when used in the Standard Model, leads to predict a proportionality 
relation, m 2 H ~ gnv R , between the squared Higgs boson mass m 2 H and the square of the 
known weak scale vr (246 GeV) through the renormalized scalar self-coupling g R ~ 1 /In A, 
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A being the ultraviolet cutoff. In this picture, the ratio m H /vR is a cutoff-dependent 
quantity that becomes smaller and smaller when A is made larger and larger. 

By accepting the validity of this picture, there are important phenomenological im- 
plications. For instance, a precise measurement of m#, say m# = 760 ± 21 GeV, would 
constrain the possible values of A to be smaller than about 2 TeV thus suggesting the 
occurrence of 'new physics' at that energy scale. 

In an alternative approach [4-6], however, this conclusion is not true. The crucial 
point is that the 'Higgs condensate' and its quantum fluctuations undergo different rescal- 
ings when changing the ultraviolet cutoff. Therefore, the relation between m H and the 
physical Vr is not the same as in perturbation theory. 

In order to remind the basic issue let us preliminarily observe that in a broken- 
symmetry phase the conventional field rescaling cannot be viewed as an 'operatorial 
statement' between bare and renormalized fields of the type, say 

" a> B (x) = Vz<f> R (x) " (i) 

In fact [7] this relation is a consistent short-hand notation in a theory allowing an asymp- 
totic Fock representation, as in QED. However, in the presence of spontaneous symmetry 
breaking it has no rigorous basis since the Fock representation exists only for the shifted 
fluctuating field. For this reason, it is the residue of the shifted-field propagator, say 
Z = Z prop , that through the Kallen-Lehmann representation is related to the normaliza- 
tion of the single-particle states. In principle, this quantity is quite unrelated to Z = Z v 
the rescaling of the vacuum field (the scalar 'condensate') which is defined through the 
physical mass and the zero-momentum susceptibility and is by no means constrained to 
be below unity. 

To be definite, let us consider a one-component scalar theory and introduce the bare 
expectation value 

VB = ($ latt) (2) 

associated with the 'lattice' field as defined at a locality scale fixed by the ultraviolet 
cutoff. Connecting to the stability analysis, such expectation value represents one of the 
absolute minima, say ip B = ±vb, of the effective potential Kff(^s) of the theory. 

Now, by Z = Z v we denote the rescaling that is needed to obtain the physical vacuum 
field 

vr = —7= • (3) 

By physical, we mean that the second derivative of the effective potential V^((pn) param- 
eterized in terms of the physical field and evaluated at <^_r = ±vr, is precisely given by 
m 2 H . Since the second derivative of the effective potential V^(v B ) is the bare zero-four- 
momentum two-point function (the inverse zero-momentum susceptibility), this standard 
definition is equivalent to define Z v as 

2 

Tn 

Z, = TT^fh = m ^ 2(0) (4) 

V eS\ V B) 

where X2(0) = 1/V^(vb) is the bare zero-momentum susceptibility. Notice that there is 
nothing in the above derivation that dictates Z v < 1. 
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On the other hand, assuming the Kallen-Lehmann representation for the shifted fluc- 
tuating field, that has a vanishing expectation value and admits a particle interpreta- 
tion, one predicts < Z prop < 1, with 'triviality' implying Z pmp — > 1 when approach- 
ing the continuum theory. Therefore, although in the standard approach one assumes 
Z v = Z pTop = 1 + O(gR), up to small perturbative corrections, one should remind the 
basically different operative definitions of the two Z's. 

For this reason, a different interpretation of triviality was proposed (see Refs. [4-6] 
and some elder works quoted therein) starting from the observations that 'triviality' does 
not require the effective potential to be a trivially quadratic function of tp. Thus, the 
theory can be 'trivial', in a technical sense, but non-trivial in a physical sense and the two 
things can coexist provided all interaction effects can be re-absorbed, in the continuum 
limit, into the vacuum structure and the physical mass of a massive free field. 

This requirement leads to consider a class of approximations to the effective potential, 
say Vcs = Kriv, where the shifted fluctuating field is governed by an effective quadratic 
hamiltonian. This includes the one-loop potential, the gaussian approximation and the 
infinite set of post-gaussian calculations where the effective potential reduces to the sum 
of a classical background energy and of the zero-point energy of a massive free field, as at 
one loop. In this class of approximations, one finds 



m H 



so that, to define the physical vr from the bare vr one has to apply a non-trivial correction. 
As a consequence, v 2 R (as now obtained from v 2 B through Eq. (J3J with a logarithmically 
divergent Z v as in Eq. (JHJ) and m 2 H ~ gRV 2 B scale uniformly in the continuum limit. At the 
same time, by construction, the shifted field becomes governed by a quadratic hamiltonian 
in the continuum limit, so that one finds Z prop — > 1 as in leading-order perturbation theory. 

By adopting this alternative interpretation of 'triviality' there are important phe- 
nomenological implications. In fact, assuming to know the value of Vr, a measurement of 
tuh would not provide any information on the magnitude of A since the ratio C = mn / Vr 
is now a cutoff-independent quantity. Moreover, in this approach, the quantity C does 
not represent the measure of any observable interaction (see the Conclusions of Ref. [7]). 

The difference between Z v and Z prop has an important physical meaning, being a 
distinctive feature of the Bose condensation phenomenon [6]. In the class of 'triviality- 
compatible' approximations to the effective potential, one finds tub/vr = 27ry / 2£, with 
< £ < 2 [6], C being a cutoff- independent number determined by the quadratic shape 
of the effective potential V e s(ipR) at <^r = 0. For instance, £ = 1 corresponds to the 
classically scale-invariant case or 'Coleman- Weinberg regime'. 

As for the standard interpretation of 'triviality', the gaussian effective potential ap- 
proach can also be extended to any number N of scalar field components [8] . In particular, 
when studying the continuum limit in the large-N limit of the theory, one has to take into 
account the non-uniformity of the two limits, cutoff A — > oo and N — > oo [9, 10]. This is 
crucial to understand the difference with respect to the standard large-N analysis. 

To check the alternative picture of Refs. [4-6] against the standard point of view, 
one can run numerical simulations of the theory and check the scaling properties of the 
squared Higgs lattice mass mf &tt against those of the inverse zero- momentum susceptibility 
1/Xiatt- According to perturbation theory, these two quantities should scale uniformly in 
the continuum limit. Therefore, a lattice computation of Z v = mf att xiatt can resolve the 
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issue. If 'triviality' is true and perturbation theory is right, assuming a single rescaling 
factor Z^ ~ Z prop , a lattice simulation has to show unambiguously that tends to unity 
when approaching the continuum limit. 

In this respect, we observe that numerical evidence for different cutoff dependencies of 
Z v and Z prop has already been reported in Refs. [11-13]. In those calculations, performed 
in the Ising limit of the one-component theory, one was fitting the lattice data for the 
connected propagator to the (lattice version of the) two-parameter form 



p 2 + m Lt 



GM = 2 T° P 2 (6) 



After computing the lattice zero-momentum susceptibility Xiatt, it was possible to compare 
the value Z v = mf att xi a tt with the fitted Z pmp , both in the symmetric and broken phases. 
While no difference was found in the symmetric phase, Z v and Z pmp were found to be 
sizeably different in the broken phase. In particular, Z prop was very slowly varying and 
steadily approaching unity from below in the continuum limit consistently with Kallen- 
Lehmann representation and 'triviality'. Z^, on the other hand, was found to rapidly 
increase above unity in the same limit. The observed trend was consistent with the 
logarithmically increasing trend predicted in Refs. [4-6]. 

Now, to our knowledge, with the exception of Refs. [11-13], there are no other system- 
atic investigations of the scaling properties of mi att vs. xutt down to values of lattice mass 
m \a,tt ~ 0.2. Therefore we might conclude that, at the present, the alternative theoretical 
scenario of Refs. [4-6] is selected by the lattice data. 

A possible objection to this conclusion is that the two-parameter form Eq. (JBJ), al- 
though providing a good description of the lattice data, neglects higher-order corrections 
to the structure of the propagator. As a consequence, one might object that the extraction 
of the various parameters is affected in an uncontrolled way thus obscuring the observed 
difference between Z^ and Z prop . 

This objection is not very serious. In fact, if 'triviality' is true, a two-parameter fit to 
the propagator data should become a better and better approximation approaching the 
continuum limit where the genuine perturbative corrections O(gn) vanish as ~ 1/lnA. 
Therefore, neglected perturbative corrections that become less and less important can 
hardly explain the observed difference between Z v and Z prop that, instead, becomes larger 
and larger. 

However, to provide additional evidence, we have decided to change strategy and 
perform a new set of lattice calculations of the zero-momentum susceptibility. In this 
way, rather than directly computing the Higgs mass on the lattice, we shall compare the 
scaling properties of xiatt with the squared mass values predicted by perturbation theory. 
Thus, at the same time, we shall be able to check: i) the previous numerical indications 
obtained in Refs. [11-13] and ii) the internal consistency of the standard interpretation 
of 'triviality' that has been accepted so far. This new computation is consistent with 
the same Z v ~ In A trend observed in Refs. [11-13] and will be presented in Sect. El 
Further evidences are presented in Sect. EJwhere a different method to combine the lattice 
observables leads to the same conclusion. Finally, Sect. 0] will contain a summary and a 
discussion of some general consequences of our results. 
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Table 1: We compare our determinations of (\<f>\) and Xlatt for given k, with corresponding 
determinations found in the literature [17]. In the algorithm column, 'S-W' stands for the 
Swendsen-Wang algorithm, while 'W stands for the Wolff algorithm. 



K 


lattice 


algorithm 


(101) 


Xlatt 


0.074 


20 3 x 24 


W 




142.21 (1.11) 


0.074 


20 3 x 24 


Ref. [18] 




142.6 (8) 


0.077 


32 4 


S-W 


0.38951(1) 


18.21(4) 


0.077 


16 4 


Ref. [17] 


0.38947(2) 


18.18(2) 


0.076 


20 4 


W 


0.30165(8) 


37.59(31) 


0.076 


20 4 


Ref. [17] 


0.30158(2) 


37.85(6) 



2 The lattice computation of 

Our numerical simulations were performed in the Ising limit that traditionally has been 
chosen as a convenient laboratory for the numerical analysis of the theory. In this limit, 
a one-component Q\ theory becomes governed by the lattice action 

Rising = -K ^ ^ X + ^M^) + ~ e M )0(x)] (7) 

x fl 

where <p(x) takes only the values ±1. Using the Swendsen-Wang and Wolff cluster algo- 
rithms we have computed the bare magnetization: 

Vb = (\4>\) , = $>Or)/L 4 (8) 

X 

(where is the average field for each lattice configuration) and the zero-momentum sus- 
ceptibility: 

Xlatt = £ 4 KM 2 MM> 2 ]- (9) 

We used different lattice sizes at each value of k to have a check of the finite-size effects. 
Statistical errors have been estimated using the jackknife. Pseudo-random numbers have 
been generated using the Ranlux algorithm [14-16] with the highest possible 'luxury'. 
As a check of the goodness of our simulations, we show in Table 1 the comparison with 
previous determinations of and xiatt obtained by other authors [17]. 

As anticipated, rather than computing the Higgs mass on the lattice as in Refs. [11- 
13], we shall use the perturbative predictions for its value and adopt the Liischer-Weisz 
scheme [21]. To this end, we shall denote by m in p Ut the value of the parameter tjir reported 
in the first column of Table 3 in Ref. [21] for any value of k (the Ising limit corresponding 
to the value of the other parameter A = 1). 

Our data for xiatt & t various k are reported in Table 2 for the range 0.1 < m input < 
0.4 (the relevant k's for m input = 0.15,0.25,0.275,0.35 have been determined through a 
numerical interpolation of the data shown in the Liischer-Weisz Table). At this point, 
we can compare the quantity 

Z V = 2/OT 4 lp utXlatt (10) 
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Table 2: The details of the lattice simulations for each k corresponding to mi nput . In the 
algorithm column, 'S-W stands for the Swendsen-Wang algorithm [19], while 'W stands 
for the Wolff algorithm [20]. 'Ksweeps' stands for sweeps multiplied by 10 3 . 



^■input 




lattice 


algorithm 


Ksweeps 


Xlatt 


Vb = (\<I>\) 


0.4 


0.0759 


32 4 


S-W 


1750 


41.714 (0.132) 


0.290301 (21) 


0.4 


0.0759 


48 4 


W 


60 


41.948 (0.927) 


0.290283 (52) 


0.35 


0.075628 


48 4 


w 


130 


58.699 (0.420) 


0.255800 (18) 


0.3 


0.0754 


32 4 


S-W 


345 


87.449 (0.758) 


0.220540 (75) 


0.3 


0.0754 


48 4 


w 


406 


87.821 (0.555) 


0.220482 (19) 


0.275 


0.075313 


48 4 


w 


53 


104.156 (1.305) 


0.204771 (40) 


0.25 


0.075231 


60 4 


w 


42 


130.798 (1.369) 


0.188119 (31) 


0.2 


0.0751 


48 4 


w 


27 


203.828 (3.058) 


0.156649 (103) 


0.2 


0.0751 


52 4 


w 


48 


201.191 (6.140) 


0.156535 (65) 


0.2 


0.0751 


60 4 


w 


7 


202.398 (8.614) 


0.156476 (15) 


0.15 


0.074968 


68 4 


w 


25 


460.199 (4.884) 


0.112611 (51) 


0.1 


0.0749 


68 4 


w 


24 


1125.444 (36.365) 


0.077358 (123) 


0.1 


0.0749 


72 4 


w 


8 


1140.880 (39.025) 


0.077515 (210) 



with the perturbative determination 

Z LW = 2kZ r (11) 

where Zr is defined in the third column of Table 3 in Ref. [21]. 

The values of Z v and Z L w are reported in Fig. 1. As one can check, the two Z's 
follow completely different trends and the discrepancy becomes larger and larger when 
approaching the continuum limit, precisely the same behavior observed in Refs. [11-13]. 
We also fitted the values for Z v to the form (A = n/a) 

Z ip = B\n{A/m R ). (12) 

Notice that the lattice data are completely consistent with the prediction Z v ~ In A of 
Refs. [4-7]. 

Of course, one might object that the discrepancy between lattice data and the per- 
turbative Zysn depends on restricting, for each given m input , to the central values of k in 
the Liischer-Weisz table. As a matter of fact, there is a theoretical uncertainty in k, for 
each given m input , that might be taken into account when comparing with the perturba- 
tive -Zlw's. For instance, for m input = 0.3, where the K-range predicted in Ref. [21] is 
0.0754(2), one might also compute xiatt for k = 0.0756 (obtaining a value xi&tt ~ 60) and 
for k = 0.0752 (obtaining a value xiatt ~ 130). As a consequence, for m input = 0.3, one 
might also conclude that Z v , as defined in Eq.flTUJ), lies in the range 0.8 — 1.7, consistently 
with the prediction of Ref. [21] Z LW = 0.932(13). 

Adopting this point of view does not represent, however, a solution of the problem. 
In fact, by inspection of Fig. 1, one can see that the range of K-values for which Z v , 
as defined in Eq.fjlUj). is still consistent with the perturbative Z L w becomes smaller and 
smaller when approaching the continuum limit. As a matter of fact, this range vanishes for 
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""imput < 0.08. This can easily be checked noticing that for m inpu t = 0.08 Ref. [21] predicts 
k = 0.07481(8), i.e. smaller than 0.0749. Therefore, by inspection of our Table 2, the 
lattice susceptibility will be definitely larger than its value for k = 0.0749, xutt ~ 1100, 
so that using Eq. fjlO)) . one gets the lower bound Z v > 1.05. This cannot be reconciled 
with the perturbative prediction, for m input = 0.08, Z LW = 0.947(10). 

In our opinion, one should not ignore this discrepancy, for instance by attempting to 
further enlarge the error bars of the perturbative predictions so as to reach a marginal 
consistency with the lattice data. In fact, in this way no meaningful test of perturbation 
theory will ever be possible. This is in contrast with the situation for the symmetric phase 
where lattice data and theoretical predictions based on the central values of k, as given in 
Table 3 of Ref. [22], agree to good accuracy, see Refs. [11-13]. As an additional check, we 
have computed the zero-momentum susceptibility for the central value k = 0.0741 that 
corresponds to m input = 0.2. From our result on a 32 4 lattice xiatt = 161.94 ± 0.67, using 
again Eq. (fTU|) . we obtain Z v = 0.960 ± 0.004 in good agreement with the corresponding 
Luscher-Weisz prediction Z LW = 0.975 ±0.010. This is also consistent with the picture of 
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Table 3: The product VgXi&tt versus k. These values correspond, for each n, to the most 
precise entries reported in Table 2. 



0.075900 3.5154 (122) 

0.075628 3.8409 (275) 

0.075400 4.2692 (270) 

0.075313 4.3674 (547) 

0.075231 4.6288 (485) 

0.075100 4.9907 (755) 

0.074968 5.8359 (622) 

0.074900 6.7349 (2186) 



Refs. [4-6] where the deviations from perturbation theory are only due to the presence of 
the scalar condensate. 



3 Further lattice determination of Z v 

There is, however, a method to analyze the lattice data where the theoretical uncertainty 
in the k — m# correlation can be eliminated. In fact, one can perform a fit of the lattice 
data with the functional form expected in the standard perturbative approach or with 
the functional form expected in the alternative scenario of Refs. [4-6]. If this is done, 
the lattice data confirm a logarithmically divergent in agreement with the previous 
indications from Fig. 1. 

To this end, let us observe that, both in perturbation theory and according to Refs. [4- 
6], the bare squared field expectation value v 2 B is predicted to diverge in units of the 
physical Higgs boson mass as 

4-~MA/m fl ). (13) 

m R 

Therefore, in the perturbative approach, where the zero-momentum susceptibility is pre- 
dicted to scale uniformly with the inverse squared Higgs mass, i.e. X2(0)m^ ~ 1, one 
expects (PT=Perturbation Theory) 

[«Ixa(0)] PT ~ln(A/m Jl ). (14) 

On the other hand, in the approach of Refs. [4-6] one predicts Z v = X2(®) m R ~ l n (A/ra^) 
so that, in this case, one would rather expect (CS=Consoli-Stevenson) 

[u|X2(0)] cs ~lB 2 (A/m fl ). (15) 

The two predictions in Eq. (jl4j) and in Eq. ffTo'j) are free of theoretical uncertainties due to 
the k — rriR correlation and can be directly compared with the lattice data for the product 
^5X2(0) reported in our Table 3. These data can be fitted to the 3-parameter form 

a|ln(K-/t c )| 7 (16) 
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Figure 2: We show the lattice data for v %Xi&tt of Table 3 together with the corresponding 
fit Eq. (JTfiJl for 7 = 2 (Consoli- Stevenson) and 7 = 1 (Perturbation Theory). 



where a is a normalization constant and we shall set the exponent 7 = 1, according to 
Eq.((T3J), or 7 = 2 according to Eq. fT5|) . 

Now, fixing 7 = 1 one obtains the totally unacceptable value chi-square~ 496 for 6 
degrees of freedom (see Fig. 2). On the other hand, fixing 7 = 2, one obtains a good 
fit of the lattice data (chi-square~ 9 for 6 degrees of freedom) with a rather precise 
determination of the critical point from the broken-symmetry phase k c = 0.074818(11) 
consistently with the estimate k c = 0.074834(15) from the symmetric phase [23]. This 
conclusion is also confirmed by the results of the full 3-parameter fit, 7 = 2. 20^21 an d 
k c = 0.07478(6) with a chi-square of 6.4 for 5 degrees of freedom. Therefore, between the 
two alternatives Eq. (fT4^) and Eq.ffTo'j). the lattice data prefer unambiguously the theoretical 
scenario of Refs. [4-6]. 

Before concluding, a brief comment about the Liischer-Weisz prediction for the product 
^5X2(0). This is based on the relation 

KX2(0)] LW = ^ (17) 
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where Zr and are given in Table 3 of Ref. [21] as functions of the mass parameter tur. 

Therefore, a comparison of Eq.(|17|) with the lattice data reported in our Table 3 re- 
introduces unavoidably the uncertainties associated with the k — itlr correlation. For 
instance, for n = 0.0759 the prediction of Ref. [21] would be [^bX2(0)] lw = 4.16 to be 
compared with the value 3.5154(122) in our Table 3. On the other hand, if one takes into 
account that the relevant values of Zr and gn in Eq. ()17|) are actually defined for the mass 
value tur = 0.4, for which the entire range is k = 0.0759(2), one might also conclude 
that the theoretical prediction lies in the range 3.76 < [v B X2(0)] LW < 4.63 that, with a 
short-hand notation, can be expressed as [v B X2(0)] LW = 4.20 ± 0.43. 

Notice, however, that the difference between the prediction 4.20 ±0.43 and the lattice 
value 3.5154(122) should not be considered a modest 1.6<j discrepancy. In fact the interval 
3.76 — 4.63 represents the entire theoretical range allowed by the numerical solution of the 
renormalization-group equations. At the same time, we observe that the agreement does 
not seem to improve approaching the continuum limit. For instance, for vtlr = 0.2, where 
the range is k = 0.0751(1), Ref. [21] predicts 5.42 < [v B X2(0)] LW < 6.31 to be compared 
with the lattice value 4.9907(755) in our Table 3. This situation might be reminiscent 
of the discrepancy in the value of Zr, pointed out by Jansen et al. [17], that was not 
improving when approaching the continuum limit. 



4 Summary and outlook 

In this paper we have reported the results of a lattice simulation dedicated to check the 
first numerical indications [11-13] for a logarithmically divergent rescaling of the scalar 
condensate. Two different methods, see Fig. 1 and Fig. 2, support the same behavior 
Z v ~ In A predicted in Refs. [4-6] from the analysis of the effective potential. 

As anticipated in the Introduction, within this scenario there is one substantial phe- 
nomenological implication that concerns the relative scaling of the physical Higgs mass 
mu and of the physical vr for which V^(vr) = m 2 H . Namely, once v R is computed from 
the bare v\ through Eq.(JHJ) with a value Z = Z v ~ B In A, so that 



2 2 1 

Vr, ~ Vr 



this v R will scale uniformly with 



R B 51nA 



™% ~ 4 (19) 



(where we have set gR ~ r4r)- 

Therefore, if vr is identified with a physical scale (e.g. vr ~ 246 GeV), there are no 
upper bounds on m# from 'triviality' since 

C = ^ ~ VAB (20) 

Vr 

is a cutoff-independent quantity [24]. The numerical value of C, however, could depend 
on the direction chosen to approach the critical line in the more general two-parameter 
(«, A) form [21] of $ 4 theory. Therefore, a new set of lattice simulations is needed to 
compute the zero-momentum susceptibility xiatt outside of the Ising limit A = 1. 
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